EXCLUSIVE GLUEBALL PRODUCTION AT THE LHC * 



M. V. T. MACHADO and M. L. L. da SILVA 

High Energy Physics Phenomenology Group, GFPAE IF-UFRGS 
Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil 
E-mail: magnus@if.ufrgs.br 



In this contribution we summarize recent results on the computation of cross sec- 
tions for glueball candidates production in quasi-real photon-photon collisions and 
on central diffraction processes, i.e. double Pomeron exchange, in heavy ion inter- 
actions at the Life are computed. In particular, we provide predictions for the 
production of exotic mesons /o(1500), /o(I710) and X(1835). The rates for these 
distinct production channels are compared. 



1 Introduction 

The gluon self-coupling in QCD opens the possibility of existing bound states of pure 
gauge fields known as glueballs (G), which are predicted by several theoretical formalisms 
and by lattice calculations [1] [2] . Many mesons have stood up as good candidates for the 
lightest glueball in the spectrum and in particular the scalar sector (J^*^ = 0++) seems 
promising. The mesons /o(1500) and the /o(1710) have been considered the principal 
candidates for the scalar glueball[3l |4]. However, in this mass region the glueball state 
will mix strongly with nearby qq states [H [S]. More recently, the BES collaboration 
observed a new resonance called A(1835) [S]. 

Recently, the clean topologies of exclusive particle production in electromagnetic in- 
teractions hadron-hadron and nucleus-nucleus collisions mediated by colorless exchanges 
such the QCD Pomeron or two photons have attracted an increasing interest [7]. The 
cross sections for these processes are smaller than the correspondent inclusive produc- 
tion channels, which it is compensated by a more favorable signal/background relation. 
Experimentally, exclusive events are identified by large rapidity gaps on both sides of the 
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produced central system and the survival of both initial state particles scattered at very 
forward angles with respect to the beam. 

In this contribution, we summarize some results presented in Ref. [8] on exclu- 
sive glueball production in two-photon and Pomeron-Pomeron interactions in coherent 
nucleus-nucleus collisions at the LHC. We also present new results for the proton-proton 
collisions at 7 TeV. In the two-photon channel, the photon flux scales as the square charge 
of the beam, Z^, and then the corresponding cross section is highly enhanced by a factor 
cx f=a 10^ for gold or lead nuclei. A competing channel, which produces similar final 
state configuration, is the central diffraction (CD) process. Such a reaction is modeled in 
general by two-Pomeron interaction. One goal of the present investigation is to compare 
the cross sections for these two channels in the production of glueball candidates. 

2 Theoretical framework and phenomenological investigations 

Let us start with the glueball candidate production in photon-photon scattering at co- 
herent heavy ion collisions using the Weizsacker - Williams approximation (EPA approx- 
imation). In such an approach, the cross section for a two quasi-real photon process to 
produce a glueball state, G, at center-of-mass energy W^^ factorises into the product of 
the elementary cross section for 77 — >■ G convoluted with the equivalent photon spectra 
from the colliding ions [7]: 



where fci_2 are the photon energies and dn/dk is the photon flux at the energy k em- 
mited by the hadron A. The photon energies determine the center-of-mass energy 
Wj.y = \/4:kik2 and the rapidity Y of the produced system. Namely, one has fci_2 = 
{Wj-y/2) exp(±y) and Y ~ (1/2) In (fci/fc2). In addition, ^/snn is the center-of-mass en- 
ergy of the ion- ion system and the Lorentz relativistic factor is given by = y/SNN/i'^rriN)- 
In the numerical calculations we use ^ysNN — 5.5 TeV and 71, — 2930 for the LHC. We 
also will present prediction for the ALICE run at ^/s = 2.76 TeV (see Table 1) 

In the EPA approximation, the flux of equivalent photons from a relativistic particle 
of charge Z is determined from the Fourier transform of its electromagnetic field. For an 
extended charge with electromagnetic form factor, Fa{Q^), the energy spectrum can be 
computed as. 




(1) 




EXCLUSIVE GLUEBALL PRODUCTION AT THE LHC 



3 



where x — k/E is the fraction of the beam energy carried by the photon and A{x) = 
1 ~ X + (l/2x^). Moreover, a = 1/137 and is the four-momentum transfer squared 
from the charge, with ~ {xmj^)'^/{l — x). 

The ghieball production in two-photon fusion can be calculated using the narrow 
resonance approximation [9^ : 

a (77 ^ G) = (2 J + 1) g r(G ^ 77) S {W^, - M^) , (3) 

where T{G 77) is the partial two-photon decay width of G, Mq is the gluebah mass 
and J is the spin of the state G. Here, we compute the production rates for the mesons 
/o(1500), /o(1710) and X(1835) [TD], respectively. The reason is due to they have been 
mentioned as possible glueball candidates by phenomenologists [TJ [2] ■ 

The predictions for the two-photon component are in practice somewhat difficult 
as the branching ratios have not been measured. To compute numerical values for the 
meson (glueball) cross section in two-photon reactions estimates for the two-photon decay 
widths are needed. The determination of them depend upon whether the meson state 
is a pure quarkonium, pure gluonic or a mixed hybrid state. A detailed discussion 
for the main fetures in each case can be found in Ref. [S], where it was verified that 
the pure glueball production gives the lowest cross section. Accordingly, the estimates 
for the two-photon width for a pure glueball resonance are r(/o(1500) — ^ 77) ~ 0.77 
eV, r(/o(1710) 77) ~ 7.03 eV and r(X(1835) ^ 77) ~ 0.021 keV. The widths 
are about three orders of magnitude smaller that for pure qq states. Therefore, as the 
two-photon cross section scales as (2 J -I- l)r(i? — >• 77), Eq. ([3]), one can consider the 
experimental feasibility of using peripheral heavy-ion collisions to determine the nature 
of the resonances discussed above. 

Let us now focus on the Pomeron-Pomeron channel. In particular, we consider the 
central diffraction (double Pomeron exchange, DPE) in nucleus- nucleus interactions. As 
a starting point we compute the DPE proton-proton cross section making use of the 
Bialas-Landshoff [Til HI] approach. It is believed that this non-perturbative approach is 
a reasonable choice due to the light mass of glueballs candidates considered and we are 
interested in the central exclusive DPE production of glueball states. In the exclusive 
DPE event the central object G is produced alone, separated from the outgoing hadrons 
by rapidity gaps, pp ^ p + gap -\-G -\- gap -I- p. In approach we are going to use, Pomeron 
exchange corresponds to the exchange of a pair of non-perturbative gluons which takes 
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place between a pair of colliding quarks. The scattering matrix is given by, 

M^Moi-j f-j Fih)Fih)expi/3{h+t2)) S^,p{V^) . (4) 

Here Mq is the amplitude in the forward scattering limit (ti = t2 = 0). The standard 
Pomeron Regge trajectory is given hy a{t) = 1 + e + a't with e « 0.08, a' = 0.25 
GeV~^. The momenta of incoming (outgoing) protons are labeled by pi and P2 {ki 
and fc2), whereas the glueball momentum is denoted by P. Thus, we can define the 
following quantities appearing in Eq. s = {pi +P2)^, si — (fci +-P)'^, S2 — (^2 + -P)^, 
ti ~ (Pi ; ^2 — {P2^k2Y ■ The nucleon form-factor is given by (t) — exp{bt) with 

b — 2 GeV^^. The phenomenological factor exp (/3 (fi +t2)) with /? = 1 GeV~^ takes 
into account the effect of the momentum transfer dependence of the non-perturbative 
gluon propagator. The factor S'gap takes the gap survival effect into account i.e. the 
probability (3^^^^) of the gaps not to be populated by secondaries produced in the soft 
rescattering. For our purpose here, we will consider S'^^p — 0.032 at y/s = 5.5 TeV in 
nucleon-nucleon collisions. 

Following [12] we find Aio for colliding hadrons, 

Mo = 32 al Dl J d^Kp^Viip^^ exp(-3 k^t^), (5) 

where k is the transverse momentum carried by each of the three gluons. V^^^ is the 
gg -> G'^ vertex depending on the polarization J of the G'^ glueball meson state. For 
the cases considered here, J — 0, one obtains the following result [T5] : 

PiVIp'2 = 7^ A (6) 

where A is expressed by the mass Mq and the width T{gg — > G) of the glueball meson 
through the relation: 

A"^ ^STTMaTigg^G). (7) 

For obtaining the two-gluon decays widths the following relation is used, F (G — )■ gg) — 
Br (G — >■ gg) Ttot{G). Here, we will be conservative and assume the resonance to be a pure 
glueball. This fact translates into an upper bound for the exclusive DPE production as 
the cross section scales with T {R gg). Following Ref. the two-gluon width can be 
computed from the resonance branching fraction in J/ip radiative decay, Jii (ip — 7G). 
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For the candidates of interest here one obtains: 

CRx\Hj(x)\''rtot Mg 

where the function Hj{{x) is determined in the non-relativistic quark model (NRQM) 
(see appendix of Ref. [M]) and cn is a numerical constant {Cr = 1, 2/3, 5/2 for 
J^^ = '", 0+"'", 2++, respectively). The masses of J/ip and of resonance are 
and Mg, respectively, and x = 1 — (Mq/M^). Based on equations above, in Ref. [M] the 
following values for the branching fractions for scalar glueballs candidates are obtained: 
Br[/o(1500)] = 0.64 ±0.11, Br[/o(1710)] = 0.52 ± 0.07. For the pseudoscalar X the 
situation is less clear due to small information on its decaying channels in radiative J/tp 
decays. Here, we set the limit case Br [X(1835)] = 1 and notice that the branching would 
be about 30 % smaller. As discussed in [8 , a consequence on the small deviation for 
the branching fraction in pure qq and glueball resonance is the difficulty in testing their 
nature using the exclusive diffractive data. 

In order to calculate the AA cross section the procedure presented in Ref. [15 is 
considered, where the central diffraction and single diffraction cross sections in nucleus- 
nucleus collisions are computed using the so-called criterion C (we quote Ref. [15 for fur- 
ther details). Using the profile function for two colHding nuclei. Tab = J (PbTA(b) Tsib— 
b), the final expression for CD cross section in A A collisions is given by jl5j : 

^SS = J d%TAA{b) exp [-A^ a^^ r^^(6)] a^^ . (9) 

where cr™ and ct^^P — S^^^ x a-p^{y/s) are the inelastic and CD cross sections in proton- 
proton case, respectively. Using Woods-Saxon nuclear densities and considering the 
inelastic cross section cr™ = 73 mb for LHC energy, ^/sAA = 5.5 TeV, we compute 
the CD cross section for nuclear collisions. The values for the inelastic cross section are 
obtained from DPMJET [16], where the scattering amplitude is parameterized using atot, 
p and elastic slope (these parameters are taken as fitted by the PHOJET model [IT]). 
We notice that for LHC energy the effective atomic number dependence is proportional 
to A^/^, which means that the nuclear CD cross section is only one order of magnitude 
larger than the nucleon-nucleon cross section. 

In what follows, we compare the two production channels and investigate the main 
theoretical uncertainties. We provide estimates of cross sections and event rates for both 
processes the LHC energies at the heavy ion and proton-proton mode. 
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3 Results and discussions 

In Table 1 numerical results for the two-photon and Pomeron-Pomeron processes are 
presented and discussed. The cross sections for glueball candidates production in photon- 
photon fusion at energy of ,/s — 5.5 TeV for lead-lead collisions are shown. The cross 
sections are sufficiently large for experimental measurement. The event rates can be 
obtained using the beam luminosity [7j: for LHC one has £pbPb = 5 • 10^^ cm~^s~^, 
which produces the following number of events. One has 3.6 • 10^, 2 • 10'^ and 4 for 
/o(1500), /o(1710) and X{1835), respectively, in the nominal LHC running time with 
ions of 10^ s (one month). We also present the results for the first ALICE run at 
2.76 TeV. For the convenience of phenomenologists we provide here a parameterization 
of a-yy{y/s — 5.5 TeV) as a function of the resonance mass. This makes simple the 
computation of event rates provided the specific meson state and its two-photon decay 
width. We obtain in the interval 400 < < 4000 MeV the parametrization: 

(Tupc [AA Rj + AA) ^ ap 
(2J+l)r(i?,;^77) 1 + {Mr/4)' 

where (Tq = 4.9147 mb/GeV and /? — —3.45335; T^j and Mn are the decay width and 
the resonance mass in units of GeV, respectively. The parameterization above allows to 
obtain estimates starting from a modeling for the two-photon width. 



Meson 




LHC (Vs = 


= 5.5 TeV) 




ALICE (Vi = 2.76 TeV) 




r^, [eV] 


Tg, [MeV] 


cr^^ [/ib] 


o-jpff [mb] 


a-y^ [/ib] 


(Tpjp [mb] 


/o(1500) 


0.77 


69.8 


1.30 


1.07 


0.78 


0.99 


/o(1710) 


7.03 


70.2 


8.60 


0.68 


4.10 


0.64 


X(1835) 


0.02 


70.3 


0.02 


0.54 


0.01 


0.50 



Table 1: Integrated cross sections for channels 77 and IP IP in A A collisions at the LHC. 

In Table 1 the results for Pomeron-Pomeron production of glueball is also presented. 
The estimates are shown for exclusive double Pomeron exchange, where estimations for 
sudakov suppression is disregarded (see discussion in Ref. [8]). For the proton-proton 
case we obtained the estimates o-jpf [/o(1500)] = 116(128)/ib, (Tff]p[/o(1710)] = 75 (83) /ib 
and o-jpff [X(1835)] = 59 (65) /ib for energy of 7 TeV (14 TeV). The Pomeron contribution 
can be separated from photon channel by imposing a cut on the impact parameter of 
collision. After imposing this kinematic cut (b > 2Ra) the Pomeron contribution is 
reduced as they are dominated by small impact parameter contributions |18) . 
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As a final remark, it is important to discuss the uncertainties on the current cal- 
culations and the experimental feasibility of detecting glueballs candidates. The main 
uncertainty here is the model dependence on obtaining the two-photon and the two-gluon 
widths for a pure glueball meson. For the two-photon width we considered a nonrela- 
tivistic gluon bound-state model of Ref. ^19,, which it could be a debatable issue and 
it is far from being optimal. There are more modern approaches as reviewed in Ref. 
[2], but this is out of the scope of present work. For the two-gluon widths, we obtained 
them from the quarkonium width based on a non relativistic boundstate calculation [T3] . 
This type of matrix elements have been discussed in Refs. [20^ giving rise to an effect of 
chiral suppression. We did not discuss the implication of those findings in present calcu- 
lation. Concerning the experimental detection, the advantage of the exclusive processes 
discussed here is clear: glueballs are probably being produced with a high cross section 
in inelastic collisions (in pp or AA reactions) but when the multiplicity is high the combi- 
natorial background is overwhelming. In exclusive production there is no combinatorial 
background. In the ultraperipheral two-photon production of glueballs, the final state 
configuration is clear: nuclei remain intact after collision and a double large rapidity 
gap between them is present (glueball is centrally produced with a low px transverse 
momenta spectrum). This type of measurement is already done at RHIC for photopro- 
duction of vector mesons and exclusive dilepton production with a signal identification 
well understood ^21,. The situation for DPE glueball production is similar, with the px 
spectrum being broader than the processes initiated by two-photons. Thus, a transverse 
momentum cut (and also impact parameter of collision) could separate the two channels. 
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